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Abstract 



The telegraph process {X(t),t > 0}, is supposed to be observed at n + 1 equidis- 
tant time points t{ = iA n , i = 0, 1, . . . ,n. The unknown value of A, the underlying rate 
of the Poisson process, is a parameter to be estimated. The asymptotic framework 
considered is the following: A n —* 0, nA n = T ^ oo as n -> oo. We show that previ- 
ously proposed moment type estimators are consistent and asymptotically normal but 
not efficient. We study further an approximated moment type estimator which is still 
not efficient but comes in explicit form. For this estimator the additional assumption 
reAjj — > is required in order to obtain asymptotic normality. Finally, we propose a 
new estimator which is consistent, asymptotically normal and asymptotically efficient 
under no additional hypotheses. 

key words: telegraph process, discretely observed process, inference for stochastic 
processes. 

MSC: primary 60K99; secondary 62M99 

1 Introduction 

The telegraph process (see Goldstein, 1951 and Kac, 1974) models a random motion with 
finite velocity and it is usually proposed as an alternative to diffusion models. The process 
describes the position of a particle moving on the real line, alternatively with constant 
velocity +v or —v. The changes of direction are governed by an homogeneous Poisson 
process with rate A > 0. The telegraph process or telegrapher's process is defined as 



where V(0) is the initial velocity taking values ±v with equal probability and independently 
of the Poisson process {N(t),t > 0}. Many authors analyzed probabilistic properties of 
the process over the years (see for example Orsingher, 1985, 1990; Pinsky, 1991; Foong 
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and Kanno, 1994; Stadje and Zacks, 2004). Di Crescenzo and Pellerey (2002) proposed the 
geometric telegraph process as a model to describe the dynamics of the price of risky assets, 
i.e. S(t) = soexp{at + aX(t)}, t > 0. where X(t) replaces the standard Brownian motion 
of the original Black-Scholes (1973) - Merton (1973) model. Conversely to the standard 
geometric Brownian motion, given that X(t) is of bounded variation, so is the geometric 
telegraph process. This seems a realistic way to model paths of assets in the financial markets. 
Mazza and Rulliere (2004) linked the process fll.lj) and the ruin processes in the context of 
risk theory. Di Masi et al (1994) proposed to model the volatility of financial markets in 
terms of the telegraph process. Ratanov (2004, 2005) proposed to model financial markets 
using a telegraph process with two intensities A± and two velocities c±. The telegraph process 
has also been used in ecology to model population dynamics (see Holmes et al, 1994) and 
the displacement of wild animals on the soil. In particular, this model is chosen because it 
preserves the property of animals to move at finite velocity and for a certain period along 
some direction (see e.g. Holmes, 1993, for an account). 

It is worth to mention that, up to now, only few references about estimation problems 
for the telegrapher's processes are known. Yao (1985) considers a the problem of state esti- 
mation of the telegrapher's process under white noise perturbation and studies performance 
of nonlinear filters. Iacus (2001) is about the estimation of the parameter 9 of the non- 
constant rate Xe(t) from continuous observations of the process. More recently, De Gregorio 
and Iacus (2006) proposed pseudo-maximum likelihood and moment based estimators for 
the telegraph process under discrete observations on a fixed time interval [0, T] when the 
process is observed with a mesh descreasing to zero. 

The aim of this paper is the estimation of the parameter A when {X(t),0 <t< T} is 
observed at equidistant times = t < . . . < t n . We assume that U = iA n , i = 0, . . . ,n, 
hence nA n = T. The asymptotic framework is the following: A n — ► and nA n = T — ► oo 
as n — > oo. 

When the telegraph process X(t) is observed continuously then N(T)/T is the optimal 
estimator of the parameter A and the statistical experiment is equivalent to the one of the 
observation of the whole Poisson process on [0,T] (see e.g. Kutoyants, 1998). This situation 
also corresponds to the limiting experiment in our asymptotic framework. 

The paper is organized as follows. Section [2] reviews some results on the telegraph process 
and presents a formula of the p-th moment of the process in explicit form. This result is 
interesting in itself because it gives new information about this model. Section [3] presents 
estimators previously introduced in the literature such as pseudo maximum likelihood es- 
timators and moment type estimators. In particular, for the moment type estimator it is 
shown that it is consistent and asymptotically normal but not efficient. Another approxi- 
mated moment type estimator is given in explicit form and it is shown that the estimator is 
consistent and (under the additional assumption nA^ — ► 0) asymptotically gaussian but still 
not efficient. Finally, Section H] presents a new estimator which is consistent, asymptotically 
gaussian and asymptotically efficient without additional assumptions. 
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2 Moments of the telegraph process 



The process X(t) is not Markovian. Conversely, the two dimensional process (X(t),V(t)), 
V(t) = V(0)(— l) N ^\ has the Markov property but a scheme of observation in which one is 
able to observe both the position and the velocity of the process at discrete time instants 
is not admissible, so statistical procedures should rely only on the observation of the X(t) 
component. The telegraph process is such that 

EX(t) = (2.1) 

and 

E * 2 « = T (< " ^) (2 ' 2 > 

(see Orsingher, 1990). Next theorem gives the general explicit derivation of the moments of 
any order of the telegraph process. In the cited reference the author mentions that derivation 
of the moments of any order can be obtained but the actual derivation was not presented. In 
some sense, next theorem completes Section 3 of Orsingher (1990) and we decided to present 
it here because it has some interest in itself. 

Theorem 2.1. Let p > 1, then 

EX(ty = (i + (-iy)(vtyf^\ 2 r (*±±\ {i^w + i^xt)} 6 — (2.3) 



Proof. We start by rewriting the p-th moment of X(t) as the sum of the two terms emerg- 
ing from the discrete and the absolute continuous part of its density (see formula (27) in 
Orsingher (1990) or (13.11) below). Therefore we have the following representation 



•rVt 

\Pp-\t 



EX{tf = J x p p(x, t)dx + (1 + (~l) p )^y — (2.4) 



-vt 



The term e xt /2 will also appear in the above integral, so we consider separately the two 
identities 

£ j +Vt x p I QvW-z 2 ^ dx = (1 + (-l) p )(vt) p (J^j 2 r (^Y^j %(At) (2.5) 

and 

f +vt j (- v / f 2 t 2 — a; 2 ) 

At / x p v y =^dx = 

^ 2.6 
2\ * fp+V 1 



-m(vt) p { (-) r(^)/^(At)-i 
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After multiplication by the factor e~ xt /2 both (12.51) and (12.61) . direct substitution in (12.41) 
gives the result of the theorem. So we need to prove the above identities and we start with 
formula ( 12.51) 

\ °° 1 / x \ 2k r+ vt 



- 2 ^(i + (-m^— j r^j %( A i ) 

= (l + (-l) P )(^) 2 (^(^JI^Af) 
where Beta(a, b) = T(a)T(b)/T(a + b), F(x) = J °° e^tf^dt, n\ = T(n + 1) and 

T ( \ 1 ( X ^ 2k+U 

In the above, the derivation of the equality 



vt , s w x ^eta (k + 1, 2+i) 

^ , „2\ K A„ _ /'I i I i\p\f„.+\p+2k+l \ 2 1 

vt 



r-p (yH 2 - x 2 ) dx = (i + {-iy){vty +2k - 



is incredibly lengthy but trivial therefore we omit it. We now calculate (12.61) 

r+vt j (A^2 t 2 _ X 2\ r+vt ~ l /A\ 2fc+1 / 9 9 9 ,fc 

A( /_„« " 7^ > dx = M L x ' £ WTw (*) " * 2 > d * 
= M f; _ 1 _ I AV t+1 (1 + ( _ in(Btr ^^(* + i.^) 

^ T(A; + 2)A;! \ ~ V 2 

,At_ /> + 1\ f / 2 \ V 1+p 



^ (1 + ( - 1)?)M YlVJ {UJ 

= (i + (-in(^|r(^) 2 /i (p _ 1} (At) - i 

because r((p + l)/2)/r((3 + p)/2) = 2/(1 □ 

Remark 2.1. From ( 12. 3D emerges that all odd-moments of the process are identically zero. 
Moreover, for x — > 0, the modified Bess el functions admit the following expansion 



x\ v I z 2 z A 



Iv ^ X ' I> + 1)U V 1 + 4(z/ + l) + 32(z/+l)(z/ + 2) + 
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from which we obtain that EX(t) p is of order t p+2 for t — ► 0. The following expansion, for 
t — > 0, will be useful in the following 

EX(t) 2 = v H 2 -lv 2 Xt 3 + ]-v 2 X 2 t 4 + o(t 4 ) (2.7) 

EX(t) 4 = vH A - -v 4 Xt 5 + -v 4 \H 6 + o(t 6 ) (2.8) 
5 5 

EX(t) 6 = 6^ 7+ 3 w + o(f8) (2 g) 

We now check that for p = 2 we recover formula (12.21) which has been derived in two 
different ways in Orsingher (1990). Indeed, for p = 2 we have 

and noticing that 

/2sinh(x) /2 x cosh(x) — sinh(x) 

iila; = W 1= — , nix) - 



direct substitution gives (12. 2p . 

Remark 2.2. The formula of the fourth moment has also a relatively simple expression, so 
we present it here. To derive the result, it is useful to know that 

12 (x 2 + 3) sinh(x) — 3xcosh(x) 
nix) = \ = 

Thus, 

3 

EX{tf = 2{vtY 2 |V5F{/| (At) + If (At)} ^ 

= 3 (j) 4 e~ A * {At (At - 3) cosh(At) + (3 + At (At - 1)) sinh(At)} 



3 Previous results on the estimation of A 

As mentioned in the Introduction, we assume that the telegraph process {X(t), < t < T}, 
with X(0) = x = 0, is observed only at discrete times < t\ < ■ ■ ■ < t n — T, with 
ti = iA n , i = 0, . ..,n hence nA n = T. We use the following notation to simplify the 
formulas: X(ti) = X(iA n ) = X«. The interest is in the estimation of the parameter A whilst 
v is assumed to be known. If the whole trajectory can be oserved, A can be estimated by 
N(T)/T where N(T) is the number of Poisson events counted in [0,T] or, the number of 
times the process switches its velocity in [0, T\. The estimation of v is always an uninteresting 
problem as, if there are no switchings in ((z — 1)A„, zA n ] then X,i — X^x = vA n , hence if 
A n is sufficiently small, there is high probability of observing N(ti+i) — N(U) = then v 
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can be estimated (actually calculated) without error. The asymptotic minimum variance 
of all the estimators for the continuous time experiment is the value of A itself because, as 
said, it is just the problem of estimating the intensity of a homogeneous Poisson process. 
We now review some estimators for this process already available in the literature and study 
some new properties of one of them. De Gregorio and Iacus (2006) considered the following 
approximated likelihood 

n 

L n {\) = L n (X\X ,X h ...,X n ) = Y[p(X i ,A n ;X^ 1 ,t^ l ) (3.1) 



i>0} 



e -XA ri 



K U n,i = 0) 



where u Uji = u n (Xi, -X»_i) = v 2 A 2 l — (Xi — Xi_i) 2 , 5 is the Dirac function and 1a is the indi- 
cator function of set A. The density p(X iy A n ; -Xf_i, appearing in (13.11) is the probability 
law of a telegraph process initially located in that reaches the position Xi at time ti. 

The above approximated likelihood is indeed the joint law of the increments Xi — which 
are considered as if they were n independent copies of the process X(A n ). The increments 
rji can be expressed as follows 

r h = X i -X i - 1 = V{0) T {-l) N ^ds = V{0){-l) N ^ T {-1) N ^- N ^ds 

and they are stationary but not independent. Conversely, the squared increments 

v f = v 2 ^ (_i)^w~^-i)d s N 2 

(or the absolute increments |^|) are independent. In their paper, the authors proposed the 
following estimator 

A n = argmaxL n (A) (3.2) 

A>0 

The estimator is proved to be unique and to exist (not so evident given the uncommon form 
of L n ) and such that A n — > N(T)/T under the condition nA n = T, A n — > as n •— >• but 
T fixed. The limiting estimator N(T)/T is the natural estimator, but X n is not consistent 
for all values of A because time T is fixed. In the same paper, the authors present numerical 
results about a least squares estimator of the following form 



In order to have consistency and hence asymptotic normality of estimators it is necessary to 
consider the asymptotics as nA n = T — > oo. We will prove that the estimator A ra is a true 
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moment type estimator which is consistent and asymptotically gaussian but not efficient. 
The estimator A n is given in implicit form and we also study an approximated moment type 
estimator which is given in explicit form and prove that it is consistent and asymptotically 
normal (under the additional condition nA\ — > 0) but still not efficient because it asymptotic 
variance is |A. A new asymptotically efficient estimator will be presented in Section HI 



3.1 The moment type estimator 

Consider the original estimator from (13.31) . The statistics 

n 

n ^— ' 

is an unbiased estimator of 



n 

8=1 



v 2 / 1 - e" 2A ° An 
Uo = /(Ao) = T [ An ~ 2A 

Observe now that /(A) is monotonic and decreasing function of A such that 

\imf(X)=v 2 Al lim/(A)=0 

A^O A^oo 

on the other side U n varies from to v 2 A 2 l , hence the minimum value of (13.31) A n is also the 
solution of 

■ " v 2 / l _ e -2AA„ 



i=l 



which means that A n is a true moment type estimator. Let Ao be the true value of the 
parameter and E and Var indicate the expected value and the variance operator under 
the true parameter A . 

Theorem 3.1. Let X n the moment type estimator solution of ( \3A\\ and let nA n — > oo, 

A n —>■ as n — > oo and let Ao > be the true value of the parameter. Then, \ n is a 
consistent estimator of Aq and such that 



^^MA n -\v)^N (o, 



where — > denotes the convergence in distribution. 

Proof. We rewrite U n as /(A n ) and study the asymptotic properties of A„ by 5-method. Let 
A = g{u) = f~ l {u) with f~ x the inverse function of /(A). Hence 

V^A^An-Ao) = \ / nA n (f' 1 (U n )-f' 1 (u )) = ^/nA n (U n -u ) — ^ - +o p ^ nA n \U n - u |) 
where 

f , m _ d .... v 2 (l + A n A + e 2 A " A (-1 + A W A)) 
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Then E y/nA n (X n — A ) = 0. Moreover, 



v 



4 



_l_l 6e 2tA (i + A n A )+e 4A " A ° (17 + 4A n A (-5 + 2A n A ))) 



Var ° ( ^ 2) " 4e^o V 



Var (\/nA n (A n - A )) 
J 1. 

A n A 2 (-l-16e 2A " A ° (l + A ri ,A ) + e 4A " A ° (17 + 4A n A (-5 + 2 A„ A ))) 



71 A " / „/, xx2 - Var (Vi ) 



4(l + A n A + e 2 ^ A ° (-l + A n A ))^ 



6A 4A 2 A ra 2 6 

~5~ + ~ + ° (A J ~> 5 A ° 



Let 



6 - ^<"? - 



in order to prove asymptotic normality, we need to prove the Lindeberg condition. Therefore, 
we need to prove that nE |£| 3 -> 0. Indeed, 



3 



3 

A/ ' / |^-/(A )| 3 + o(A 2 )) 



%W(A )) J 

3 

AqAh , I 
+ o(A; 



□ 



3.2 An approximated (but explicit) moment type estimator 

Consider again (12.21) . Some algebra, or Remark 12.11 give the following expansion 



V 2 / 1 - p- 2AA « 

E (X, - X,_0 2 = T- A, 



A V " 2A 



v^f 2XA n - l(-2XA n f) - |(-2AA n ) 3 ) + o(A 3 ) 
A l An 2A 



= v 2 Al - ^ 2 AA 3 + o(A 3 ) 
Therefore, an approximated moment type estimator is the following 



and A* is a weighted sum of the independent random variables rjf. Remark that v 2 A 2 n — 
(Xi — Xj_i) 2 is exactly zero if no Poisson event occurs in the time interval (ti, U + i\. This fact 
will be used to evaluate expected values of related quantities in the following. This estimator 
is qualitatively not different from the estimator X n in equation (13.31) but the cost of having 
an estimator in explicit form is paid by the need of use of the additional condition nA\ — ► 
in order to obtain asymptotic normality. This hypothesis has been also used in the "high 
frequency" sampling for discretely observed diffusion processes (see e.g. Florens-Zmirou, 
1989; Yoshida, 1992). 

Theorem 3.2. Let Ao be the true value of the parameter. Then, under the condition A n — ► 0, 
nA n — > oo as n — > oo the statistics A* in (I3.5P is consistent estimator of X . Moreover, under 
the condition nA^ —>■ it is also asymptotically Gaussian, i. e. 

V^(K-x )^n(o^x 

Proof. Consistency is trivial, indeed 

E q A; = Ao + o(l) 

We now study its asymptotic distribution, in particular we are interested in the following 
random variable 

y q n ( n ■) n 

V^~n(K - Ao) = £ V 'K ~ ( X i - - o AoA 3 „ \ = £ & 

We now show that 

e {v^a;(a;-Ao)} ~ v^Af 

holds true, where a ~ b means "a is the same order of 6". This is where the assumption 
nA^ — > is needed in order to obtain asymptotic normality for this estimator. 



< n\d\ (P{N((ti-i,ti)) = °} + PW(*i-i,*;]) = 1} + PiNttU-uti]) > 1}) 



3n 



2^AI 
nA^ 



- 2 -X Ay + A A„ UaI - - ^X Al ) + o(Af 



'nAl -> 



Moreover, Cov(£j,£j) = 0, for i ^ j. Now we need to prove Lindeberg condition. We 
calculate the variance of £j making use of the following two expansions derived from (12. 7p . 
(BHD and (EH]) 

Er? 2 = E(X, - X^f = v*Al - 2 -XAl + \v^Ai + (Al) 
E^ 4 = «^ - ^ 4 AA* + ^A 2 A* + o{Al) 
Vvt = v*Al - 6 -v«XAl + ^ 6 A 2 A* + o(Al) 
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so that 

Var(tf) = E oV f - (E ^) 2 = ^X Al - |, 4 A^Af t + o(A f 6 J 



hence 



Thus 

Var { v^A^A; - A )} = nVar(&) = ^A - ^AqA„ + ^o(A n ) - ^A (3.6) 
Now, let S n = 71 + 7 2 + • • • + 7n and = VarS^ with 

3 I .,2*2 /v v \2 2 \ A3 I Si 



7. = ^TT v'K ~ (* - X^Y - -A A^ 



2ra; 2 A 3 [ n v 1 l ~ L > 3 u n J 
We have that 

2 w o ^t- / \ 1 6 A 8A 2 , 9 o(A r 



4 = VarS n = nVar( 7l ) = -^Var(&) 



A n ~" 5nA„ 5 n 4t> 4 nA„ 
and we observe that 

Eo{l(N>«-}7?}= / lndPo= f 7^dP <7^ / N 3 dP < ^ E oN 

{|7i|>"n} {|7i|>"n} {\~ii\>£Sn} 



thus 



1 v ^ , 21 ra 

— 2^ E {l{\ lz \>es n }li ) < — 

Now we maj orate E 



n i=i n 

13 



E|7t|3 = 8^3aI E 

27 



VA^-j^A*) - ^ 2 



8v 6 n 3 A? t 



2 

3 AoA " 



PjiV^.!,^) = 0}+ \KA 2 n + o(A 2 n )\ 3 P{N[t^ 1 ,t i ) > 1} 



^8JAl(^ + A„|^ + o(AD| 3 ) 



1 1 

^ -7 + 



n 3 n 3 A 2 

in the above we have majorated P{N[ti-i,U) = 0} by 1 and denoted by K a generic constant. 
Easy manipulation shows that 

4 = 4\/^ = — - — 
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hence 

nE | 7 ,| 3 < AI + ^ ^ 

with the notation a < b indicating "a is majorated by something of the same order of 6". 
Hence we have proven Lindeberg condition and asymptotic normality is established. Trivially 
and using exactly the same arguments, the Lindeberg condition can also be proved directly 



4 An asymptotically efficient estimator 

In the previous section we have seen that the estimators A* and A n are not efficient because 
their asymptotic variance is |A which is strictly greater than A which is the asymptotic 
variance of N(T)/T, the asymptotically efficient estimator of the continuosly observed Pois- 
son process. By a different approach we now present an asymptotically efficient estimator. 
Consider the following statistic 

Y n 1 n 

~ Xn = Yl 1 {lwK«A n } = — r- Yl 1 {-V([*i- 1 A))>1} ( 41 ) 

ILL\ n . "-'-in • , 

1=1 1=1 

The statistic A n is not a good estimator of A for fixed A n . Indeed, 

I ™ i _ e -A A„ 

E (A„) = — — 2jE (l{|^| <tJ A„}) 



i=i 



Hence, we propose the following estimator 



A n = --^log(l-A n A n ) (4.2) 



A.. 

and next theorem proves that it is the efficient estimator in this context. 

Theorem 4.1. Let A n — ► 0, nA n —>■ oo as n —>■ oo. Then, for all fixed Ao > 0, the estimator 
An in (14.21) is consistent, asymptotically normal and attains the minimum variance, i.e. it 
is asymptotically efficient 

V^A^(A n -A ) ^iV(0,A ) 

Proof. In order to prove consistency and asymptotic normality of A n we first prove the same 
properties for A n . We have seen that A n in (14. ip is such that 

I — g-AoA„ 

E (A„) = = A + A^A* + o{A 2 n ) - A 
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hence consistency of A n follows trivially. We now prove asymptotic normality Let us consider 
the following quantity 

U n = \/ nA n (\ n — E A n 
1 n 

= /T77\ — Yl { 1 {lwl<«A n } - E o (l{|lJ,|<t>An}) } 
V il'^n i=1 

1 ™ 



1=1 



i=l 



with 



6 = ^{l te |<^}-(l-e- AoA ")} 



We have that E ^ = thus E (£4) = 0. Moreover, 

Var (&) = Vax (l{| % |<„ A „}) 

= E (l{| 7?i |<„A„})(l - E (l{| % | <uA „})) 
__ H Aq A„ Aq A n 

= A A n + o(A A n ) 

hence 

Varo(C^) = — — n(A A n + o(A A n )) = A + o(l) 
nA„ 

Finally, the £j's are independent because they only involves the absolute value of the incre- 
ments t]i, hence once we prove the Lindeberg condition we also have asymptotic normality. 
For large n it holds true that < l/y/nA n , then Lindeberg condition is trivially true 

n 

EMl{|6l> e >£}->0 

i=i 

and the following result holds true 

U r AN(0,\ ) 

Now we need to prove asymptotic normality of A n in (14.21) . Let 

/(«) = log(l - «A n ), /'(u) = = - J—, 

A n dii 1 — «A„ 

and 

/~ 1 — e - ^ ^ 

An = /(A n ), Aq = / Aq = T 



12 



then, by the so-called 5-method, 

V^K(X n - A ) = y^A~ n (f(X n ) - /(A )) 

= \/nA n (\ n - A )/'(Ao) + o v {\fnA n \\ n - A |) 

= \fnE^(\ n - A )- 7— r— + o p (l) 

1 — A A n 

hence 

V^lAn-Ao) ^iV(0,A ) 

□ 

5 Conclusion 

In this paper we have shown that a consistent and asymptotically efficient estimators of the 
parameter A of the telegraph process observed at discrete times can be obtained under the 
high frequency sampling when the length of observation interval [0, T] increases. We have 
also shown that the moment type estimators A* and A n are consistent but not efficient. They 
are also asymptotically Gaussian (the approximated moment type estimator A* requires the 
additional assumption nA^ — > 0). This paper also contains the explicit formula for the 
moments EX(t) p of any order p > 1 of the telegraph process and a simple expansion for 
small t. 
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